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Very high energy gamma-rays are expected to be absorbed by the extragalactic background light 
over cosmological distances via the process of electron-positron pair production. Recent observations 
of cosmologically distant gamma-ray emitters by ground based gamma-ray telescopes have, however, 
revealed a surprising degree of transparency of the universe to very high energy photons. One 
possible mechanism to explain this observation is the oscillation between photons and axion-like- 
particles (ALPs). Here we explore this possibility further, focusing on photon- ALP conversion in 
the magnetic helds in and around gamma-ray sources and in the magnetic held of the Milky Way, 
where some fraction of the ALP flux is converted back into photons. We show that this mechanism 
can be efficient in allowed regions of the ALP parameter space, as well as in typical configurations of 
the Galactic Magnetic Field. As case examples, we consider the spectrum observed from two HESS 
sources: lESllOl-232 at redshift z = 0.186 and H 2356-309 at 2 = 0.165. We also discuss features 
of this scenario which could be used to distinguish it from standard or other exotic models. 

PACS numbers: 95.85.Pw, 98.70.Vc 98.70.Rz 14.80.Mz 
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I. INTRODUCTION 

Very high energy (VHE, > 0.1 TeV) gamma-rays 
are expected to be attenuated over cosmological distances 
by infrared, optical and ultraviolet photons of the ex- 
tragalactic background light (EBL) via the process of 
pair production (7vhe + 7ebl — ^ e+e"). This mecha- 
nism is expected to strongly suppress the VHE spectra 
from distant sources and opens the possibility of mea- 
suring the spectrum and density of the EBL via gamma- 
ray observations Recent findings by imaging atmo- 
spheric Cherenkov telescopes, however, indicate a rel- 
atively large degree of transparency of the universe to 
gamma-rays 0, Hj. Taken at face value, the data seem 
to require a lower density of the EBL than expected 
and/or considerably harder injection spectra than ini- 
tially thought 0, 1^. On the other hand, it is also pos- 
sible that some exotic mechanism is responsible for the 
observed lack of attenuation. 

In Ref. 6|, it was argued that the oscillations of pho- 
tons into axion-like particles (ALPs) in extragalactic 
magnetic fields could provide a way of avoiding the expo- 
nential suppression of the gamma-ray spectrum resulting 
from pair production. This mechanism is especially ef- 
fective for extremely distant sources and in the presence 
of relatively strong (nanogauss scale) extragalactic mag- 
netic fields. While consistent with present bounds, there 
is no observational evidence for such strong fields, which 
appear hard to generate dynamically. For example, re- 
cent simulations based on magneto-hydrodynamic ampli- 
fication of seed magnetic fields driven by structure for- 
mation predict fractions of extragalactic space contain- 
ing nanogauss or stronger fields ranging from less than 
0.1% to 10% Also, recent results from the Pierre 
Auger Observatory appear to favor magnetic fields which 



are closer to the weaker estimates, with upper limits for 
Mpc coherence lengths in the 0.1 1 nG range in most of 
the sky 0,[l3]- In this respect, although photon- ALP os- 
cillations in extragalactic magnetic fields may in principle 
be possible, it would require unsupported assumptions on 
the fields in order to act as an effective mechanism for 
avoiding the attenuation of VHE gamma-ray sources. 

A somewhat different approach was taken in Refs. [ll| 
and [l^l where it was argued that the magnetic fields 
present in astrophysical accelerators themselves would 
lead to significant photon- ALP conversion over a large 
range of ALP parameter space, including a small re- 
gion of the very interesting parameter space in which 
Peccei-Quinn axions provide a solution to the strong 
CP problem. This effect is expected to become sig- 
nificant at gamma-ray energies, including the range 
being explored by the ground based gamma-ray tele- 
scopes HESS ilAMAGIC [13, VERITAS [l|] and 
CANGAROO-III[ll|, as weh as of the forthcoming sateL 
lite mission, GLAST Tt^. In this article, we explore the 
possibility of reducing the cosmological attenuation of 
VHE gamma-rays through a photon- ALP mixing mech- 
anism, but invoking only known magnetic fields: those 
needed to confine and accelerate cosmic rays at the 
gamma-ray sources, where a significant ALP flux is pro- 
duced via oscillations. To reconvert a fraction of the ALP 
fiux back into photons, well after the primary photon flux 
has disappeared due to the absorption on the EBL, we 
rely only on the magnetic field of the Milky Way. This 
mechanism can lead to the appearance of a substantial 
fiux of gamma-rays (up to ^30% of the low-energy ex- 
trapolation) at energies where none is expected to be 
seen. It has also some other very peculiar observational 
predictions, which we will also discuss. 

The remainder of this paper is structured as follows. In 
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Sec.|TT]we describe the mechanism proposed. In particu- 
lar, we devote Sec. Ill Al to the details of the reconversion 
of ALPs in the Galactic Magnetic Field. In Sec. IIIBI we 
describe how wc treat the absorption of photons by the 
EBL during the propagation from their sources to the 
Milky Way. In Sec. Ill C|, in order to illustrate the impact 
of such mechanism in a realistic scenario, we present the 
reconstructed spectrum for two HESS sources: lESllOl- 
232 at redshift z = 0.186 and H 2356-309 at z = 0.165. 
In Sec. lIIII we discuss and summarize our results, focusing 
on the distinctive observational features of the scenario 
outlined here. 



II. THE MECHANISM 



As discussed in Refs. [ll|, [l^, given the typical sizes 
and magnetic field strengths present in astrophysical ac- 
celerators, significant photon to ALP conversion can oc- 
cur in or near the VHE gamma-ray sources over a large 
range of allowed ALP parameter space. In the limit 
of complete "depolarization" of the photon-ALP sys- 
tem, one expects 1/3 of the original photons to be con- 
verted into ALPs at the source above a critical energy 
E = m1^/{2ga-yB), where rua is the ALP mass, ga-y is 
the ALP-photon coupling and B the field strength. That 
is, for £' ^ 5 there could potentially be an ALP flux 
from VHE gamma-ray sources as large as ^50% of the 
residual gamma-ray flux exiting the source. The energy 
£ naturally falls in the gamma-ray band. Although it is 
not crucial for the viability of the mechanism being dis- 
cussed, we shall assume that 5 <C 0.1 TeV, thus insuring 
that the spectrum observed by ground based gamma- 
ray telescopes will not demonstrate peculiar spectral fea- 
tures resulting from the onset of photon-ALP oscillations. 
As we shall see below, this condition is naturally ful- 
filled over the most viable range of p arameters for which 
the mechanism is efficient (see [ll[ for details). Thus, 
over the whole range of energy observed by atmospheric 
Cherenkov telescopes, the ALP flux generated at the 
source via oscillation follows exactly the same spectrum 
of the photons at the source. 

This photon-ALP mix then propagates towards the 
Milky Way without further oscillations, since in absence 
of an external field the ALPs and the photons effec- 
tively decouple. During propagation over cosmological 
distances, the spectrum of photons is depleted via pair 
production. As a result, upon reaching the Milky Way, 
the flux will be dominated by ALPs, some of which can 
be converted back into photons by the Galactic Mag- 
netic Field. Since this step is crucial and it constitutes 
the more original part of the mechanism we propose, in 
Sec. Ill A] we describe it in detail. 



A. Reconversion in the Milky Way 

In a given direction in Galactic coordinates, 
(6,/), the unit vector, s, with components 
{cos&cos/, — cos6sin/, sin&} can be constructed. 
Through the usual orthogonalization procedure of 
Gram-Schmidt, § can be completed into a K"^ basis 
(along with, say, a and b, two vectors orthogonal to s). 
One can thus decompose the magnetic field, B(x), into 
components {Ba^ Bjj, Bs\ in this basis and show that 
the probability of an ALP converting into a photon (of 
any polarization) while traveling a distance s along s is 
given by p^ : 



/'ds'e'(^"-^p')^'B„(s') 

JO 



(1) 



This probability is a factor of 2 larger than P^_,a(s) 
given in Ref. [l^ due to the two photon states avail- 
able for the ALP to convert into. In this expression, the 
quantities entering the phase factor are — —rn^/2E 
and Api — uj'^y/2E, where Wpi = Aira ne/we is the 



plasma frequency, E the photon energy, rrie the electron 
mass, Ue the electron density, and a the fine structure 
constant. For the following considerations, it is use- 
ful to introduce the following dimensionless quantities: 
.gii = 507/10^" GeV"\ Skpc = s/kpc, TOncV = m/neV 
and Ejw = E/TeV. Results from the CAST experi- 
ment [1^ provide a direct bound on the ALP-photon cou- 
pling of gu < 8.8 for rua < 0.02 eV, nominally below the 
long-standing globular cluster limit [20|. For ultra- light 
ALPs {ma < 10~^^ eV), the absence of gamma-rays from 



SN 1987A yields a stringent limit of gu < 1 [2l| or even 
gii < 0.3 [22]. In the range 10-"eV < Wa < lO'^eV, 
the CAST bound is the most general and stringent. For a 
range of /leV-scale masses, bounds from ADMX (2^ are 
stronger, but rely on the assumption that axions consti- 
tute the Galactic dark matter. In these illustrative units, 
the phases can be rewritten as: 



AaSkpc = -0.77 X 10- 



'ncV 



TcV ' 



Apiskpc = -1.11 X 10-7 {^)s^pcE 



(2) 



ToV ■ 



The integrals in Eq. ([T|) are very small unless the phases 
of the integrands vanish over typical Galactic sizes of 
~ 10 kpc. Barring fine-tuning, we must require |Ai x 
lOkpcj <C 1, leading to 



< VI. 3 X IO^-BtoV, 



Tie < 0.90 X 10^ cm-3 EtcV 



(3) 



For the energies of interests for atmospheric Cherenkov 
telescopes, the second condition is always satisfied in the 
Milky Way. The first condition only holds if nia ^ 
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10~^ eV. The optimal range of parameters for this mecha- 
nism is thus 10~^'^eV < nia < 10"^ eV, over which cou- 
pUngs as large as gn ~ 8 are consistent with present 
bounds. 



The conversion probability also depends on the geome- 
try of the Galactic Magnetic Field, which is not very well 
known, especially in the directions toward the Galactic 
Center and away from the Galactic Plane. In Fig. [l] we 
show isocontours in ALP-photon conversion probability 
for three different models of the Galactic Magnetic Field 
(see Ref. [13 )• Following the authors' initials, we label 
the models as PS TT [26|, and HMR Phys- 
ically, the PS model includes a dipolar field which pro- 
duces a small component perpendicular to the Galactic 
Plane in the Solar neighborhood but is responsible for a 
strong field in the Galactic Bulge. The PS model also 
includes a toroidal halo field at kpc distances above and 
below the Galactic Plane. The HMR model does not in- 
clude a dipolar field, but has a more prominent field in 
the halo characterized by the same symmetry pattern of 
the field in the disk. Finally, the TT model is the most 
conservative of the three, with a weaker intensity for the 
local coherent field and no prominent regular field in the 
Galactic Center or in the halo. Each of the models imple- 
ment a thin disk field following the spiral pattern of the 
Galaxy, but the symmetry with respect to the Galactic 
Plane in the TT model is opposite to the one assumed in 
HMR and PS. Further details can be found in Ref. [H. 
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The isocontours of ALP-photon conversion probability 
shown in Fig. [1] were calculated using Eq. for gu = 5 
and assuming that the conditions of Eq. ^ are satisfied. 
This means that the probability can be rescaled accord- 
ing to gfi, as long as Pa^^ <C 1 (Eq. ^ indeed only 
holds at leading order in perturbation theory). In the 
limit of full mixing in the Galaxy, one would obtain a 
reconversion probability of 2/3. 

For illustration, we also show in Fig. [T] the positions 
of all currently known VHE gamma-ray sources at a red- 
shift of 0.1 or greater, symbol-coded according to Table I. 
It is clear that for gn — 5 there are regions of the sky in 
which the reconversion probability can be 20% or even 
larger, and that a significant fraction of the sky corre- 
sponds to a probability larger than 10%. It is intriguing 
to notice that, in the HMR model, many of the VHE 
gamma-ray sources lie within or nearby these regions. We 
conclude that appreciable reconversion probabilities are 
possible, although difficult to predict given the scarcity 
of our knowledge regarding the structure of the Galactic 
Magnetic Field. In the remaining of this paper, we shall 
assume reconversion probabilities of ~0.1 and explore the 
resulting phenomenological consequences. Note that for 
gii ^ 5, nia ~ 10~^eV and microgauss-scale fields, the 
critical energy £ falls in the sub-GeV range, consistently 
with the assumption put forth in Sec. II. 
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FIG. 1: The probability of ALP-photon conversion in three 
different models of the Galactic Magnetic Field (PS, TT, 
HMR from top to bottom), for gn = 5 and assuming that 
the conditions of Eq. ([Sjl are satisfied. Also shown are the 
locations of the known very high energy gamma-ray sources 
at redshift greater than 0.1, labeled according to the symbols 
found in table I. 
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Object 


z 


I 


h 


Ref. 


Symbol 


3C279 


0.536 


305.10 


57.06 


[28] 


A 


PG 1553+113 


> 0.25 


21.91 


43.96 


[29] 


▲ 


lESlOll+496 


0.212 


165.53 


52.71 


m 


□ 


lES 0347-121 


0.188 


201.93 


-45.71 


m 


• 


lESllOl-232 


0.186 


273.19 


33.08 


[32] 


■ 


1ES1218+304 


0.182 


186.36 


82.73 


[33] 


+ 


H 2356-309 


0.165 


12.84 


-78.04 


[34] 


* 


lES 0229+200 


0.140 


152.94 


-36.61 


[35] 





H 1426+428 


0.129 


77.49 


64.90 


[36] 


o 


PKS 2155-304 


0.117 


17.73 


-52.25 


[3J 


X 



TABLE I: The current catalog of known very high energy 
gamma-ray sources at redshift greater than 0.1. The locations 
of these objects are also shown in Fig. [1] 



B. Photon absorption onto the EBL 




E [TeV] 



FIG. 2: Suppression of the gamma-ray spectrum due to pair 
production for sources at redshifts of z = 0.1, 0.2, 0.3, 0.4, 
and 0.5, from top to bottom. 



To study the quantitative predictions of our model 
and make comparisons with the available data, we must 
also account for the standard dimming of the VHE pho- 
tons from cosmologically distant sources. The spectra 
of distant VHE gamma-ray sources will be attenuated 
through electron-positron pair production. This attenu- 
ation takes the form exp(— r), where t(E^ zq) is given by 
the equation: 



l-z / , dn , , _ , ^ 

Aijj — {uj,z)cj{E,uj,z), (4) 



where 



(l + z)i?(z)y dw 



-(1 - ii)a^^ 



{E,uj,n)dn. (5) 



Here, Zq is the source redshift, H{z) is the rate of Hubble 
expansion, E and u are the (appropriately redshifted) 
source and background photon energies, respectively, fi 
is the cosine of the angle between the incoming and target 
photon, and cr-^^ is the cross section for electron-positron 
pair production [38j . 

The distribution of EBL photons in energy and red- 
shift which we use, dn/doj, is based on the integrated 
galaxy light model of Ref. [s^ and the data of Ref. [i^ . 
We assume that the shape and normalization of the back- 
ground radiation does not change except for redshifting 
between the source and the observer, which as suggested 
by the models of Ref. [llj is a reasonable assumption for 
z < 0.4. For the sources being considered here (see ta- 
ble I), we expect this to be a reasonable assumption. In 
Fig. [21 we plot the attenuation as a function of observed 
energy for sources at various redshifts. Note that it is 
not only the overall intensity of the source but the shape 
of the spectrum that is affected by this suppression. 



C. Results 

At this point, we have all of the ingredients needed to 
compute the spectra of VHE gamma-ray sources, includ- 
ing ALP reconversion. In Fig. [3l we show the results for 
a source spectrum of dA^/d£' cx E~^ with F = 2 after 
propagating from z = 0.2 (top) and z = 0.5 (bottom). 
For comparison, we also include the results for a stan- 
dard scenario with no ALPs. Clearly, in the presence of 
ALP-photon mixing, the quasi-exponential cutoff is inter- 
rupted, and a plateau in the spectrum appears at high 
energies. 

We next use this model to compute the source spectra 
of two HESS sources, lESllOl-232 at z = 0.186 and H 
2356-309 at z — 0.165, using the data given in Ref. [2|. 
Our results are shown as Fig. H) The points shown with 
error bars represent the spectrum as measured by HESS, 
whereas the other points denote the source spectrum (be- 
fore propagation) required in order to obtain the observed 
spectrum. From these figures, it is clear that extremely 
hard injection spectra are required (F « 0.6 and -0.2 
for the two objects, respectively) to match the observed 
spectrum after propagation if one restricts oneself to the 
fiducial astrophysical models for the EBL. If the effects 
of ALP-photon oscillations are included, however, the 
required slope of the source spectrum is softened to a 
reasonable value for each source (F « 2.0 and 1.8, re- 
spectively in this example, assuming Pa-y = 0.1 in our 
Galaxy). In each frame, we also show the case in which 
the EBL is reduced to only 45% of its fiducial value. 
This is the scenario proposed by the HESS collaboration 
to explain the observed lack of attenuation. 

With this exercise, we have shown that the effects of 
ALP-photon mixing in known magnetic fields enables one 
to fit the VHE data with a reasonable spectral index at 
the source (F w 2). This mechanism can accommodate 
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FIG. 3; The gamma-ray spectrum from a source with an in- 
jected spectrum of dN/AE oc , after propagation over a 
distance of 2 = 0.2 (top) and z = 0.5 (bottom). Results are 
shown with and without the effect of photon- ALP oscillations. 
The ALP-photon mixing mitigates the impact of absorption 
via pair production and leads to a plateau in the spectrum at 
high energies. We have calculated the effects of ALP-photon 
mixing assuming a source conversion probability of 0.3 and a 
Milky Way reconversion probability of 0.1. The vertical axis 
is in arbitrary units. 



FIG. 4: The points shown with error bars represent the 
observed gamma-ray spectrum of H 2356-309 (top) and 
lESllOl-232 (bottom), from the data of Ref. 2]. The other 
points shown denote the spectra required to be injected at 
the source in order to generate the observations after ac- 
counting for the absorption by the EBL. Using a conven- 
tional EBL spectrum and density, the reconstructed source 
spectra are required to be extremely hard {AN/AE oc 
and AN/AE oc E"'^). The mechanism of photon- ALP mixing, 
however, can soften the required spectral slope at injection to 
acceptable values. See text for more details. 



the typical predictions for the EBL without need to mod- 
ify the infrared to ultraviolet emission models from cos- 
mic star formation, allowing spectra at the source that 
are certainly much easier to obtain in typical acceler- 
ation models. Also, this model has several distinctive 
phenomenological consequences, which we shall comment 
upon in the next section. 



III. DISCUSSION AND CONCLUSION 

In recent years, imaging atmospheric Cherenkov tele- 
scopes have discovered many new TeV gamma-ray 
sources, some of which are cosmologically distant (see 
Table I). These sources provide us with a useful probe 



of the extragalactic background light in the infrared to 
ultraviolet range. In particular, the spectra of very high 
energy gamma-ray sources at cosmological distances are 
expected to be attenuated by this background through 
the process of electron-positron pair production. Re- 
cent observations seem to indicate a far greater degree of 
transparency of the universe to very high energy gamma- 
rays than previously estimated [1, [3| . While astrophys- 
ical explanations of these observations have been dis- 
cussed 0, Q, it has also been proposed that this lack 
of suppression could be the result of photon mixing with 
axion-like particles (ALPs), with oscillations occurring in 
the presence of the magnetized intergalactic medium [6,] . 

In this article, we have discussed an alternative mech- 
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anism involving photon- ALP mixing, but requiring only 
known magnetic fields: namely the fields in or near the 
gamma-ray sources which are needed to confine and ac- 
celerate particles (which in turn are necessary to pro- 
duce gamma-rays), and the Galactic Magnetic Field 
of the Milky Way. We have shown that this mecha- 
nism can be efficient for an ALP mass in the range 
lO^^^eV <ma< 10"^ eV and couplings of 511 ~ 4. A 
precise prediction of the modification expected for a given 
source is precluded by our current ignorance regarding 
the detailed structure of the large scale magnetic field of 
the Milky Way. A very robust prediction of this model, 
however, (which is very different from both the standard 
expectation and the model proposed in Ref. 0) is that 
the degree of dimming observed is expected to be depen- 
dent on the galactic coordinates {l,b}. Thus, the most 
striking signature would be a reconstructed EBL density 
which appears to vary over different directions of the sky: 
consistent with standard expectations in some regions, 
but inconsistent in others. 

The absorption of gamma-rays via pair production 
does not affect only point-like source emission, but also 
the diffuse (or unresolved) spectrum (see e.g. Ref. [43|). 
Thus another signature of this mechanism might be an 
anomalous and direction-dependent spectral behavior of 
diffuse radiation. This will be challenging to detect with 
present atmospheric Chercnkov telescopes, however, as 
they are not well suited for the study of diffuse radia- 
tion over large fields-of-view. In contrast, the satellite- 
based experiment GLAST will observe the entire sky, al- 
though with much smaller exposures than are possible 
with ground based gamma-ray telescopes. Provided that 
the characteristic energy for the onset of the ALP-photon 
conversion mechanism naturally falls in the energy range 
to be explored by GLAST, the study of the diffuse radi- 
ation offers an independent test of the crucial ingredient 
of this model, namely the role of the Galactic Magnetic 
Field. 



Through the conversion of photons into ALPs in the 
Galactic Magnetic Field, a peculiar, direction-dependent 
depletion of photons should arise in the Galaxy (at up to 
the 30% level). This idea was considered previously for 
a different range of the ALP parameter space, impacting 
the diffuse X-ray background [1^. The anisotropics in 
the diffuse emission should anti-correlate with the regions 
of stronger fields as detected, for example, in Faraday ro- 
tation maps. Put another way, a peculiar dimming of the 
diffuse radiation in the 1-100 GeV range may be detected 
from the same regions where stronger than expected TeV 
sources are present. 

A detailed test of such a scenario will become increas- 
ingly possible as a more complete picture of the high 
energy gamma-ray sky is developed. Since existing at- 
mospheric Cherenkov telescopes have a very small field 
of view, surveys of large fractions of the sky are unfeasi- 
ble. A targeted survey may be possible, however, as new 
sources in the GeV range are discovered by GLAST. Ul- 
timately, next generation gamma-ray observatories, such 
as CTA [il] or AGIS [11] , will enable considerably more 
detailed studies. 

The mechanism described here provides yet another 
example of the new opportunities for discovery made pos- 
sible as a result of the recent progress in the field of high 
energy gamma-ray astrophysics. 
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